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Abstract 

The diversity product and the diversity sum are two very important parameters 
for a good-performing unitary space time constellation. A basic question is what the 
maximal diversity product (or sum) is. In this paper we are going to derive general 
upper bounds on the diversity sum and the diversity product for unitary constellations 
of any dimension n and any size m using packing techniques on the compact Lie group 
U{n). 

1 Introduction 

Let A be a matrix with complex entries. A* denotes the conjugate transpose of A. Let 
denote the Frobenius norm of a matrix, i.e., 



\\A\\ = y/tl (AA*). 

A square matrix A is called unitary if A* A = AA* = I, where I denotes the identity matrix. 
We denote by U(n) the set of all n x n unitary matrices. U(n) is a real algebraic variety and 
a smooth manifold of real dimension n 2 . For the purpose of this paper a unitary space time 
constellation (or code) V is simply a finite subset of U (n) , 

V = {A 1 ,A 2 ,-- - ,A m }cU(n). 

We say V has dimension n and size m. Unitary space time codes have been intensely studied 
in recent years and we refer the interested readers to E3 El and the references of these 
papers. The readers will find the motivation and engineering applications of such kind of 
codes. The quality of a unitary space time code is governed by two important parameters, 
the diversity product and the diversity sum. 



*Both authors were supported in part by NSF grants DMS-00-72383 and CCR-02-05310. The first author 
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Definition 1.1. The diversity product [5] of a unitary space time code V is defined through 



min{| det(A - B)\»\A, B <E V, A ^ B}. 



The diversity sum ^T] is defined as 



5> 



1 



min{\\A- B\\\A,B e V,A^ B}. 
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V is called fully diverse if JJV > 0. As explained in a space time code with large 
diversity sum tends to perform well at low signal to noise ratios whereas a code with a large 
diversity product tends to perform well at high signal to noise ratios. A major coding design 
problem is the construction of unitary space time codes where the diversity sum (or product) 
is optimal or near optimal inside the set of all the space time codes with the same parameters 
n, m. We would like to remark that for every positive integer n and m, a Haar distributed 
random space time code is fully diverse with probability 1. A simple proof can be found 



The purpose of this paper is to derive for n and m tight upper bounds for the diversity 
product Y\V and the diversity sum J^V. When n = 1 then trivially | det(A — B)\ = \\A — B\\ 
and it follows that V = n ^ i n this situation. The following lemma states that for every 
space time code V, V is an upper bound for Yl V and by having an upper bound for Yl V 
we immediately also have an upper bound for Yl V. The readers can find the statements 
about the relationship between Yl V and V in [TT] , for completeness we include a detailed 
proof. 

Lemma 1.2. For any unitary space time code V, 



where U, V are unitary matrices and Cj > for j = 1, 2, • ■ ■ , n are the singular values of C. 
First we are going to prove 



If Cj = for some j, then the inequality is trivial. Hence we assume Cj > for all j's. 
Because U, V are unitary matrices, it follows that 



in0. 




Proof. Let C be an n x n complex matrix with singular value decomposition 



C = [/diag(ci,c 2 , • • • ,c n )V, 



dfit(C# < -=||C|| 




Similarly one verifies that 
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Applying Cauchy-Schwarz inequality, we have 



n 



En . 



n 



-|det(C)|» < — ^, 
2 1 v ;| - 2Jn' 



Hence one concludes that for an n x n square matrix C, 

-\\c\l 

By the definition of f| V, ^ V and the above inequality one gets 



□ 



Of course it would be desirable to know for every n and m what the largest possible value 
of ^2 V is. This is the motivation of the following definition. 

Definition 1.3. Let A(n,m) be the infimum of all numbers such that for every unitary 
space time code V of dimension n and size m, one has 



^2 V < A(n,m). 



Remark 1.4. As pointed out by Liang and Xia ^T] there exists a constellation V of dimen- 
sion n and size m with V = A(n, m). This is due to the fact that U(n) m is a compact 
manifold. 

The exact values of A(n, m) are only known in very few special cases. In the case n — 1, 

one checks that A(l, m) = sin — for m > 2. When n > 2 and m — 3, one has A(n, 3) = 
When m = 2, we have A(n, 2) = 1 for n > 2. For n = 2, the following values were computed 

in mi- 
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Liang and Xia [TT] observed the connection between a unitary constellation and an Eu- 
clidean sphere code and beautifully derived an upper bound for 2 dimensional unitary con- 
stellations which is very tight when m < 100. In this paper we present a new general upper 
bound for A(n, m) for every dimension n and every size m while improving certain results 
in [TI]. To the best of our knowledge the new upper bounds we derived are tighter than any 
previously published bounds as soon as m is sufficiently large. 



2 Upper Bound Analysis 

In this section we are going to study the packing problem on U(n) and derive three up- 
per bounds for the numbers A(n, m). All the resulted bounds are derived by differential 
geometric means and all bounds can be viewed as certain sphere packing bounds. 
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From a differential geometry point of view we can view U (n) as a n 2 -dimensional compact 
Lie group. U(n) is also naturally a submanifold of the Euclidean space M? n . In this way 
U (n) will have the induced geometry of the standard Euclidean geometry of M? n . Finally 
there is a third way to see U(n) as a submanifold of another Riemannian manifold S(n) and 
we will say more later. 

The basic strategy for computing the upper bounds for A(n, m) is as follows. Given a 
unitary space time code V = {A\,A2,--- ,A m }, around each matrix Aj we can choose a 
neighborhood N r (Aj) with radius r (the radius will be specified later). Let Vj = V(N r (Aj)) 
be the volume of the neighborhood N r (Aj). If all the neighborhoods are non-overlapping, 
then necessarily we will have 

m 
3=1 

where V(U(n)) denotes the total volume of unitary group U(n). This inequality in turn will 
result in an upper bound for the numbers A(n, m). By employing different metrics (Euclidean 
or Riemannian) and by considering different embeddings of U(n), we derive three different 
upper bounds for A(n,m). 

Let Mi be the manifold consisting of all the n x n Hermitian matrices, i.e. 

Mi = {H\H = H*}. 

Mi has dimension n 2 and can be viewed isometrically as Euclidean space R™ 2 . Assume that 
H = (Hjk) and assume that Hj k = Xj k + iyjk- We use (dH) to denote the volume element 
of A^i, where 

(• \ n(n— 1)/2 n n 
% - J l\ dH u /\ dH jk /\ dH ]k = /\ dx u /\ dx jk /\ dy jk . (2.1) 

' 1=1 j<k j<k 1=1 j<k j<k 

With a small abuse of the notation, one can check that the volume element of M2, the 
manifold consisting of all the n x n skew-Hermitian matrices, can be written as 

(. \n(n— 1)/2 / 1 \ n n n 
\) (-) /\dH u J\dH jk /\dH jk = l\d Vll /\dx jk /\dy jk . (2.2) 

' ^ ' 1=1 j<k j<k 1=1 j<k j<k 

For a unitary matrix U, if we differentiate U*U = I, we will have 

U*dU + dU*U = 0. 

Therefore U*dU is skew-Hermitian. The following lemma will characterize the volume el- 
ement of U(n). For the terminologies in this lemma, we refer to a standard differential 
geometry or integral geometry book, e.g. [81 112j. 

Lemma 2.1. The volume element ofU(n) induced by the Euclidean space R 2 ™ 2 is bi-invariant 
and the volume element can be written as (U*dU) up to a scalar constant. 

Proof. The bi-invariance comes from the orthonormality of U(n). (U*dU) is left-invariant 
according to the definition. Indeed for a fixed yet arbitrary unitary matrix V, 

(VU)*d(VU) = U*V*VdU = U*dU. 

Since U(n) is a compact Lie group and any compact Lie group is unimodular, (U*dU) is also 
right-invariant. Because the bi-invariant n 2 differential forms are unique up to a scalar, one 
concludes that the volume element can be written as (U*dU). □ 
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The following theorem will represent the volume element of U(n) in another way. One 
will see that it is closely related to the eigenvalues of unitary matrices. 

Theorem 2.2. For the Schur decomposition of a unitary matrix 0: 

6 = f/diag(e*\e ie V- - ,e i6n )U*, (2.3) 

we will have 

{Q*dQ) = Y[ \e ie > - e idk \ 2 d6 1 A d0 2 A ■ ■ • A d6 n A (U*dU - diag(CTW)). (2.4) 

j<k 

Proof. Let D = diag (e 101 , e td2 , ■ ■ ■ , e t6n ) and take the differential of Equation (|2.3|) . 

dQ = dUDU* + UdDU* + UDdU*. 

It follows that, 

0*dO = UD*U*dUDU* + UD*dDU* + UdU* = U{D*U*dUD + D*dD)U* + UdU*. 

Due to the right-invariance of the volume element in U(n), it follows that 

(Q*dQ) = (U*Q*d<dU) = (D*U*dUD - U*dU + zdiag (dO u d6 2 , • • • , d6 n )). 

Note that (D*U*dUD - U*dU) jk = (e^ - e Wk )U jk , therefore the diagonal elements of 
D*U*dUD — U*dU are all zeros and the off diagonal elements are scaled version of the 
ones of U*dU. According to formula ([2.2)1 . the claim in the theorem follows. □ 

The following theorem calculates the volume of a small neighborhood with Euclidean 
distance r. Because of the homogeneity of U(n), the center of this small "ball" is chosen to 
be / without loss of generality. For a unitary matrix U, we assume e j, s are its eigenvalues, 
i.e., U ~ diag (e l6>1 , e 2 , • • • , e t0n ). For a fixed unitary matrix A, let 

Uj?(n,A) = {U eU{n)\\\U - A\\ < r}. 

Again because of the homogeneity of U(n), V(Uf (n, A)) does not depend on the choice of 
A. In the sequel V(U^(n)) will be used to denote V(Uf(n, A)) for any unitary matrix A. 
Let S(n) denote a 2n 2 — 1 dimensional sphere centered at the origin with radius y/n, i.e., 

S(n) = {(x 1 ,x 2 ,--- ,^2n 2 )ki + ^H \-x 2 2n 2 =n}. 

Apparently U(n) is a submanifold of S(n). For a particular point So G S(n), let 

S r (n,S ) = {Se S(n)\\\S - S \\ < r}. 

Theorem 2.3. Let 

Di = {(9u 9 2 , ■ ■ ■ , n )\ - 7T < B 5 < it for j = 1, 2, ■ • • , n} (2.5) 



and 



then 




D,= {(B 1 ,B 2 ,-.. ,^)|>; s in 2 |<^[>, (2.6) 



fL nn n^fc - e^\ 2 dB x de 2 ■■■d9 n 
V(U?(n)) = JJ T D l! lj< *\ _ jai2 \ n X - - V(U(n)). (2.7) 



JJ Dl U,< k \e^-e^\ z d9 1 d0 2 ---d6. 
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Proof. Note that \\I — U\\ 2 < r is equivalent to ^^^ =1 sin 2 < For a given unitary 
matrix 6, the Schur decomposition = U* di&g (e 1 ® 1 , e 102 , ■ ■ ■ ,e ie ")U is unique if OjS are 
strictly ordered. So if we take the integral of formula ()2.4j) over the integration region 
disregarding the order of Oj's, we will obtain n\ times the volume of V(Uf(n)). Thus the 
volume of U^(n) will be 



V(U r E (n)) = 




(U*dU -diag (U*dU)). 



Using the same argument, we will derive the volume of U(n): 

V{U{n)) = i- II TT \e ie * - e Wk fdO^ ■ ■ ■ d9 n 1 1 {U*dU - diag (U*dU)). 

rijJD.f^ ' JJu{n) 

Compare the two derived volume formula, the claim in the theorem follows. □ 
Remark 2.4. By the Weyl denominator formula j3] one can replace 

II TT | e ^ - e i6k fde x de 2 ■■■de n 

1 j<k 

with (2n) n n\. We keep it as it is to make the formula literally understandable. 

There are several approaches to derive upper bounds for the diversity sum. The first 
approach considers U(n) as a submanifold of S(n), then chooses the non-overlapping neigh- 
borhoods to be small balls with radius r (with regard to the Euclidean distance). This will 
result in the first upper bound (Bl) which we derive in this paper. 

Theorem 2.5. Let D\ and D 2 be defined as in (|2.5|) and ()2.6|) . Assume r§ = 7q (n, m) is 

the solution to the following equation (with variable r): 

mil FT |e^ - e i9k fd6 1 d6 2 ■ ■ ■ d9 n = II TT |e^ - e i0k \ 2 d6 1 d6 2 ■ ■ ■ d6 n , (2.8) 

JJD 1 nD 2j<k JJ Dx j<k 

then 

I ( r E\2 (rE\4 

A(n,m)<^^-^J_. (Bl) 

Proof. For a fixed yet arbitrary unitary constellation V = {Ai,A 2 ,--- ,A m }, consider m 
small non-overlapping neighborhoods S r (n, Aj) in S(n). We can increase r such that there 
exist I, k such that S r (n, A{) and S r (n, AfA are tangent to each other. Apparently 

UffaAj) = Sri^A^nUin), 

for any j. Since SV(n, A,-)'s are non-overlapping, we conclude that U^(n, A?)'s are non- 
overlapping. Therefore we have 

m 

^F(C/ r E (n,A,))<F(C/(n)), 

3=1 

6 



that is 

mV(Uf(n)) < V(U(n)). 

One can check that V(Uf(n)) is an increasing function of r, so any r satisfying the above 
inequality will be less than the solution to the equality: 

mV(U?(n)) = V(U(n)), 

which is essentially Equality ( |2.8|) . So we conclude that r < rjf. 

Note that any two points So, Si G S(n) with two non-overlapping neighborhoods S r (n, S ) 
and S r (n, Si) will have distance \\Sq — S±\\ > — r 4 / (An), where the equality holds only 

if S r (n, So) and S r (n, Si) are tangent to each other. Apply the argument to A/s and note 
that Ai and Ak are the closest pair of points with \\Ai — Ak\\ = 2yV 2 — r 4 /(4n), we reach 
the conclusion of the theorem. □ 

For a fixed So G S(n), consider S r (n, Sq) C S(n). Let r = r(n, r) denote the maximal 
number r such that S r (n, Si), S r (n, S 2 ), ■ ■ ■ , S r (n, S T ) are non-overlapping and S r (n, Sj) is 
tangent to S r (n, S ) for j = 1,2, ■■• ,n. One checks that r(n,r) does not depend on the 
choice of Sq. In this sense r(n,r) can be viewed as generalized kissing number j3] on an 
Euclidean sphere. For a fixed n dimensional unitary constellation V = {Ai, A%, ■ ■ ■ ,A m }, 
let r(V) denote the maximal radius r such that S r (n, A x ), S r {n, A 2 ), ■ ■ ■ , S r (n, A m ) are non- 
overlapping. Let r opt = r opt (n,m) denote the maximal r(V) over all possible n dimensional 
unitary constellation V with cardinality m. One checks A(n, m) = r opt (n,m) / \/2n. The 
following theorem and corollary give a lower bound for the optimal diversity sum A(n,m). 

Theorem 2.6. Let Di be defined as in ()2.5|) and assume that Tq = r^n, m) is the solution 
to the equation (12. 8J) . Let 




< (r(n, r opt (n, m)) + 1) // J] - | 2 ^^ 2 ■ • ■ d9 n . 

J JD^Dz j<k 

Proof. According to the derivation of r$, we have 

mil | \ \e ie > - e i6k \ 2 d6xd6 2 ■ ■ ■ d9 n = 1 1 \\ |e^ - e i6k fd9id9 2 ■ ■ ■ d9 n . (2.9) 
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Assume that V = {Ai, A2, • • • , A m } is an n dimensional unitary constellation reaching 
r opt (n,m), i.e., r(V) = r opt (n,m). For simplicity let r = r(V). Let m! denote the maximal 
number such that S r (n, Ai), S r (n, A 2 ), • ■ ■ , S r (n, A m ), • • • ,S r (n,A m /) are non-overlapping. 
Let 7"i = 2\Jr 2 — r 2 / (An), we claim that 

m' 

U(n)c\JuZ(n,Aj). 

3=1 

Otherwise suppose there is a unitary matrix A £ Uj=i Uf^n, Aj), then ||A — > r\ 
(see Theorem 12. 5|) . Thus S r (n,A ) does not intersect with S r (n,Aj) for j = 1,2, ••• ,m'. 
Therefore one can find m' + 1 small balls with radius r which are non-overlapping. This 
contradicts the maximality of ml. Thus we have YlT=i ^(^ri( n J A?)) — ^(^( n ))' that is 



J| |e^ - e idk \ de l( ie 2 ■■■d6 n > j \ |e^ - e id *\ dQ x dQ 2 ■ ■ ■ d6 n . (2.10) 
We further claim that 

m! < (m- l)(r(n,r) + 1). (2.11) 
By contradiction assume that m' > (m — l)(r(n, r) + 1) + 1. Let 

tang(j) = {/|1 < / < rn',S r (n,Ai) tangent to S^n, Aj)}. 

According to the definition of r(n, r), we know the cardinality of tang (j) is less than r(n, r). 
We first pick ji from {0, 1, • • • , m'}, then pick j 2 from {0, 1, • • • , m'} — tang (ji). And we 
continue this process by always picking j^+i from 

k 

{0,1,- •• ,m'}-|Jtang(iO- 

Since the cardinality of the above set is strictly greater than when k < m — 1, we can pick 
jii hi ' ' ' j jm from the index set {1, 2, • • • , m'} such that S r (n, AjA, S r (n, Aj 2 ), • • • , S r (n, Aj m ) 
are non-overlapping and every two of them are not tangent to each other. Then we can find 
a small enough real number e > and increase the radius r to r + e such that 

S r +e( n , Aj-j, S r+£ (n, Aj 2 ), ■ • • , S r+£ (n, Aj m ) 

are still non-overlapping. However this contradicts the maximality of r = r opt (n,m). 
The combination of the three formulas ()2.9|) . ()2.10|) . (|2.11|) will lead to 



ff Dl U ]<k \e^-e^\ 2 de 1 d9 2 ---d9 n 



Jf Dl nB 3 n ]<k \e^-e^\ 2 d6 1 d6 2 ...d6 n 

( Un FLv* \e i6j - e l ^fdQ x dQ 2 ■■■dQ n 
< JJ ^ LL ^<k\ LJ_J — -1 I Mn,r) + 1). (2.12) 

~ V IJ Dl nn 2 U j<k \& ~ e^\ 2 de ld 9 2 ■ ■ ■ d0 n ' 

Note that the inequality above is in fact stronger than the claim in the theorem. We can 
reach the conclusion of the theorem by relaxing the right hand side of the inequality (by 
ignoring —1). □ 



S 



Corollary 2.7. When m — > oo ; asymptotically we have 

A(n,m) > 2^nr E (n,m)]-(T{2n 2 - 1) + l)~ 1/n2 . 

Proof. We only sketch the idea of the proof. Intuitively U E (n,Ao) looks more "flat" when 
m — > oo (consequently r — > 0), so V(U E (n, A)) can be approximated by the volume of 
U E {n, A)'s projection to the tangent space of U(n) at A : 

TT | e ^ _ e iek \ 2 d9 1 d9 2 ■ ■ ■ d9 n ~ C(rf )™ 2 
for some constant C. The same argument will lead to 



// 



J J \ e *i _ e ^ | 2 ^^ 2 • • • d9 n ~ C(2r opt )" 2 

for the same constant C. For any fixed n, r(n, r) will approach to the standard kissing num- 
ber in Euclidean space r{2n 2 — 1) when r goes to zero. Combining the three approximations, 
we reach the claim according to the previous theorem. 

□ 

U(n) is a compact Lie group equipped with a Riemannian metric. Given two points 
A , Ai G U(n), one can always find a geodesic j(t) (mapping from [0, 1] to U(n)) which will 
connect these two points, i.e. 7(0) = Aq and 7(1) = A\. Recall that the Euclidean distance 
of Aq and A\ is defined to be ||A — We further define the Riemannian distance between 
Aq and Ai to be: 

dist (A ,A 1 ) = [ h\t)\\dt. 
Jo 

As a Lie group U (n) is homogeneous. In particular one has that 

dist (A), A^ = dist (UAq, UA x ) = dist (A U, A X U) 

for any U G U(n). The following theorem utilizes the homogeneity and the relationship 
between the Riemannian distance and Euclidean distance to derive another upper bound for 
the diversity sum in general and it is the base of the second approach. 

Theorem 2.8. Let /(•) and g(-) be two fixed monotone increasing real functions. If 

g(\\A - AH) < dist (Aq, A,) < f(\\A - All) 

for any two unitary matrices A and A\, then 

A(n,m)<g- 1 (2f(rZ(n,m)))/(2V^). 

Proof. For a fixed unitary constellation V = {A, A, • • • , A m }, consider 

U r E (n, A,), U r E (n, A), • • • , U r E (n, A m ) 
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for r > 0. We can increase r until there exist j and k such that Uf(n,Aj) and Uf(n,Ak) 
are tangent to each other at a point A . As examined in Theorem 12.51 one can make a 
conclusion that r < r$ '(n,m). Accordingly we have 



dist (Aj, A k ) < dist (A j7 A ) + dist (A k , A ) 

</(||A i -Ao||) + /(|K-A || 

On the other hand since (7 is monotonically increasing one has: 



2f(r)<2f(r*(n,m)). 



\\Aj -A k \\<g ' (dist (Aj,A k )). 
The combination of the above two inequalities will lead to 

\\Aj-AkW < g- 1 (2f(r^(n,m))). 

Immediately we will have 

J2v<g-\2f(r*(n, m )))/(2^). 

Since V is an arbitrary unitary constellation, the claim in the theorem follows. □ 

Based on the above theorem, the following corollary gives the second upper bound (B2). 

Corollary 2.9. For a real number r, let \r\ denote the greatest integer less than or equal 
to r, then 



A(n, m) < sin 



\ 



7T- 

n 



r o) 2 ( n i m ) 



H — arcsin 

n 



r o) 2 ( n i m ) 



{tq ) 2 (n, m) 



:b2) 



Proof. Consider / and another point U = Vdiag (e ldl , e %d2 



)V*, where -tt < dj < 



it. 



It is known that 0j the geodesic from I to U can be parameterized by 



7(t) = Vdiag (e l0l \ e 



iOtt je 2 t 



„i9 n t 



)V* 



where < t < 1. The Riemannian distance from / to U is 



dist (i,u) = y/ei + q + --- + e*. 

We want to derive <?(•)>/(') as i n Theorem 12.81 Suppose the Euclidean distance between / 



• 2 , • 2 ^2 . . . 2 @n 2 I A 
sin h sin !-••• + sin — = r /4. 



and U is r, i.e., 

2 ' 2 ' ' 2 

After substituting with Xj = sin 2 Oj/2 and denoting = arcsin 2 y/x, we convert the above 
problem to the following optimization problem: 

Find the minimum and maximum of the function 



F( Xl ,x 2 , ■ ■ ■ , x n ) = e\ + e 2 2 + ■ ■ ■ + e 2 n = a(g( Xi ) + G(x 2 ) + ■■■ + G{x n )) 
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with the constraints x\ + x 2 + • • • + x n = r 2 /4 and < x 3 < 1 for j — 1, 2, • • • , n. Since G(x) 
is a convex function on [0, 1], we derive the lower bound of F(x±,x 2 , • • • , x n ), 

4narcsin 2 (r/(2v / ^)) < F(x±, x 2 , ■ • • ,x n ). (2.13) 

In the sequel we are going to calculate the upper bound of F(xi,x 2 , ■ ■ ■ ,x n ). Without loss 
of generality, we assume < x\ < x 2 < • • • < x n < 1. Let k = [r 2 /4j and a = r 2 /4 — k, we 
claim that F(x±,x 2 , • • • , x n ) will reach its maximum when 

{0 1 < j < n - k - 1 
a j = n — k 
1 n- k + l< j <n 

Suppose by contradiction that F reaches its maximum at (x±,x 2 , • • • , x n ) with X\ > 0. Now 
from 

xi + x n _ k + x n _ fc+ i H h x n < r 2 /A = k + a, 

surely one can find x' n _ k , x' n _ k+1 , ■ ■ ■ ,x' n such that 

X\ -\- X n —k -\- X n _fc_|_x -|- • • • + X n X n _j ( -\- X n _j t _^i -\- ■ ■ ■ -\- x n , 

with x'j > Xj for j = n — k, n — k + 1, ■ ■ ■ , n. Now set x\ = 0, x* = Xj for j = 2, 3, • • • , n — k — 1 
and x* = x', for j = n — k, n — k + 1, • • • , n. By the mean value theorem, there exist (/s 
with x\ = < (1 < xi and Xj < Q < x* for j — 2, 3, • • • , n such that 

n 

F(x{,x 2 ,--- ,x* n ) -F(xi,x 2 ,--- ,x n ) = ^2G'(Q)(x* - xj). 

Since G(x) is a strictly convex function, we have 

0<G / (Ci)<G / (C 2 )<---<G'(Cn). 

Now 

n 

F(xl,x* 2 , •••,<)- F( Xl ,x 2 , ■■■,x n )> G"(C 2 )(5>* - xj)) - G"(Ci)(x! - x\) 

3=2 

= (G"(C 2 ) - G'(Ci))(^i - xj) = (G'(C 2 ) - G'(Ci))x! > 0. 

This contradicts the maximality of F at (x±, x 2 , • • • , x n ). Applying exactly the same analysis 
to x 2 , x 3 , ■ ■ ■ , x n _fc-i, x n ^k we deduce that Xj = for j = 2, 3, • • • , n — k — 1 and = a. 
So the upper bound of F can be given as 



F(x ± ,x 2 , ■ ■ ■ , 




Take g(r) = 2y / narcsin(r/(2y / n)) and f(r) = 2a/A;7t 2 /4 + arcsin 2 ^fa, the corollary follows 
according to the previous theorem. □ 



11 



0.5 i r 1^5 2 2 4 6 r 8 10 12 14 

Legend Legend 

****** upper bound (B1) (n=3) upper bound (B1) (n=1 00) 

upper bound (B2) (n=3) upper bound (B2) (n=1 00) 



Figure 1: The comparisons of two upper bounds as functions for n = 3 and n = 100 



Note that both upper bound (Bl) and upper bound (B2) depend on r§(n, m). In Figured] 
we plot both upper bounds as functions of r§(n,m) for 3 and 100 dimensions. One can see 
that if and only if tf(3, m) > 2.0881, the upper bound (B2) is tighter than the upper bound 
(Bl). While for the 100 dimension case, the upper bound (Bl) is tighter than the upper 
bound (B2) if and only if r® (100, m) > 11.9155. In fact it can be checked that asymptotically 
when n is large enough, upper bound (B2) is tighter than upper bound (Bl) if and only if 
tf(n,m) > 1.1892^. 

For a packing problem on a manifold, alternatively one can choose the neighborhood to 
be a small "ball" with Riemannian radius r. This will be our third approach to derive an 
upper bound for the diversity sum. For a particular A 6 U(n), let 

U?(n,A) = {U e I7(7i)|dist (U, A) < r}. 

Note that the constraint dist (U, I) < r is equivalent to 

e{ + el + • • • + e 2 n < r 2 . 

Therefore we apply the same argument as in the proof of Theorem 12.31 and conclude that: 

V(U r R (n)) = j ^! nD4llj<fcl J — -V(U(n)), 

II Dl U J<k \e^-e^\ 2 d9 1 d9 2 ...d9 n 

where Di was defined in (12 .5|) and 

n 

D 4 -.= {(e 1 ,e 2 ,---,e n )\Y,^<r 2 }. (2.14) 

3=1 



m ) 1 



Instead of considering the Euclidean neighborhoods Uf(n, Ax), U^(n, A 2 ), ■ ■ • ,U^(n,A, 
we can consider the Riemannian neighborhood U^(n, Ai), U^(n, A 2 ), • • • , U^(n, A m ). Utiliz- 
ing the fact that the Euclidean distance ||A-; — Aj^\ and the Riemannian distance dist (Aj, A^ 
are related (compare with Formula ()2.13|0 : 

4n arcsin (\\Aj - A k \\/(2y/n)) < dist (A h A k ) 



12 



for any two unitary matrices Aj and Ak, we can derive the third upper bound (B3). The 
proof of the following theorem is very similar to the one of Theorem 12.81 and for the sake of 
brevity we omit it. 

Theorem 2.10. Let D\ and be defined as in (12. 5|) and f!2.14|) and assume r^(n, m) is 
the solution to the following equation (with variable r): 



m 



FT |e^ - e i6k \ d9 x d9 2 ■ ■ ■ d6 n = FT \e*' - e i6k \ dQ x dd 2 ■ ■ ■ d6 n , (2.15) 



then 



r^(n. in ; 



n 



A(w,m) < sin u ' . (B3) 



We gave three approaches to derive upper bounds for the diversity sum and hence also for 
the diversity product. All of them involve the calculation of rjf (n, m) or j-q (n, m), which are 
the solutions of equation (|2.8|) and equation (|2.15|) . respectively. Fortunately we are dealing 
with finding a root of a monotone increasing function (recall that both V(Uf(n, m)) and 
V(U^(n,m)) are monotone increasing functions with respect to r), the bisection method |2] 
will be highly effective to solve this kind of problem. Our numerical experiments for small 
size constellations with small dimensions show that upper bound (B3) is looser than the first 
two upper bounds. However when m goes to infinity, these three upper bounds give almost 
the same estimation. This makes sense because asymptotically the small balls look like a 
n 2 dimensional ball in Euclidean space. One can see the derived upper bounds for 2 and 3 
dimensional constellations in Figure |21 

We compare the derived upper bounds with the currently existing one presented in (llj . 
For n = 2 the upper bounds derived by Liang and Xia JT] tend to be better when m < 100 
and our bounds become tighter when m > 100 (see the following Table Q). For n > 3 Liang 
and Xia jTTj outlined a method by considering a sphere packing computation in S(n). It is 
our belief that this method will result in a weaker bound than the upper bounds we derived 
in this paper. For the sample programs to do the upper bound calculation, we refer to 



Table 1. For n = 2 the following table compares the upper bounds in JT] with our new 
bounds (Bl) and (B2). 



m 


24 


48 


64 


80 


100 


120 


128 


1000 


upper bounds in [11 J 


0.6746 


0.6193 


0.5969 


0.5799 


0.5632 


0.5499 


0.5452 




upper bound (Bl) 


0.7598 


0.6603 


0.6131 


0.5932 


0.5578 


0.5425 


0.5347 


0.3270 


upper bound (B2) 


0.7794 


0.6734 


0.6235 


0.6026 


0.5654 


0.5496 


0.5415 


0.3285 



One interesting fact about the limiting behavior of A(n, m) (when m — > oo) is its connec- 
tion to the Kepler problem 0. Certainly one can use Kepler density to obtain a tighter 
bound of the diversity sum asymptotically. 



3 Conclusions and Future Work 

We presented three approaches to derive upper bounds for the diversity sum of unitary 
constellations of any dimension n and any size m. The derived bounds seem to improve the 
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Figure 2: Upper bounds for 2 and 3 dimensional constellations 

existing bounds when n = 2 and m > 100. When n is large the exact computation of r$ 
is rather involved and hence it is also computationally difficult to compute the bounds (Bl) 
and (B2). Nonetheless it is our belief that the resulting upper bounds (Bl) and (B2) become 
fairly tight as soon as m is sufficiently large. 

It was pointed out that the resulted upper bounds also apply for the diversity product, 
although the bounds seem to be less tight in this situation. The future work may involve the 
derivation of a tighter upper bound analysis for the diversity product of unitary constellations 
using differential geometric means. 
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